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1. Introduction 

Digital in-line holography (DIH) is widely used in microscopy for biological applications [IJE], in 
3D Holographic Particle Image Velocimetry (HPIV) for fluid mechanics studies [3] and in refrac- 
tometry [4]. In most theoretical DIH studies, the optical systems or the objects used, for example 
particles, are considered to be centered on the optical axis 5,6 . Nevertheless, in many practical ap- 
plications, the systems and objects are not necessarily centered. This is not a problem under plane 
wave illumination. But more and more studies involve a diverging beam, or astigmatic beams directly 



obtained at the output of fibered laser diodes or astigmatic laser diodes. The position of the object 
in the beam becomes particularly relevant. To reconstruct the image of an object, a reconstruction 
parameter must be determined in order to find the best focus plane. For wavelet transformation, 
the parameter is the scale factor [7]. For Fresnel transformation, the parameter is the distance z 
from the object to the quadratic sensor plane [5] and for the fractional Fourier transformation, the 
parameters are the fractional orders [9lll0| . It is considered here that the same parameter can be used 
for all positions of the object. Recently, an analytical solution of scalar diffraction of an elliptical 
and astigmatic Gaussian beam (EAGB) by a centered opaque disk under Fresnel approximation has 
been proposed. By using the fractional Fourier transformation, a good particle image reconstruction 
is obtained However, there are no recent publications in DIH theory demonstrating that recon- 
struction parameter is constant for all transverse positions of the object in the field of the beam. 
This is, however, particularly important if we want to develop wide-angle metrologies. 

In this publication, the aim is to demonstrate that the same fractional orders can be used to 
reconstruct an image of the particle whatever its position in the field of the beam. We exhibit the 
effect of the Gaussian beam on the reconstructed image. In the first part of this publication, the 
model of the analytical solution of scalar diffraction of an EAGB by a centered opaque disk is 
revisited to take into account a decentered object. In the second part, the definition of the fractional 
Fourier transformation is recalled and this transformation is used to reconstruct the image of the 
particle. It is in this part that we demonstrate that the same orders can be used. Finally, we propose 
to illustrate our results by simulations and experimental results. 

2. In-line Holography with an elliptic and astigmatic Gaussian beam 

The basic idea in DIH is to record by a CCD camera the intensity distribution of the diffraction 
pattern of an object illuminated by a continuous or pulsed wave [12] . Figure |T]) represents the 
numerical and experimental set-up where all parameters are identified. The incident Gaussian beam 
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of diameter ui crosses a plano-convex cylindrical lens. This cylindrical lens acts only in the maxis. Its 
focal length is equal to f v = 200mm. After propagation in free space, the elliptical and astigmatic 
Gaussian beam illuminates an opaque particle. In the plane of the object, the beam widths along 
the £-axis and 77-axis are defined by and w^. The wavefront curvatures of the astigmatic beam are 
denoted by R q with q = ^,rj. The CCD camera is located at a distance z from the opaque particle. 

The basic model to describe the intensity distribution of the diffracted beam by an object recorded 
by the CCD camera is the integral of Kirchhoff-Frcsncl given by the scalar integral for the complex 
amplitude A in the quadratic sensor plane: 

A = eM [^ Z) I £ T (^)exp(f [^- x f + {r] -yf]) dt dri , (1) 

l\z J R 2 \\Z V 

in which -Et(£ 5 v) ' 1S the product of the optical incident beam, denoted E(£, 77), by the spatial trans- 
mittance of the shifted opaque 2D-object, denoted 1 — — Co, V~ Vo)- The quadratic sensor records 
the intensity defined by \A\ 2 . If we consider that the function Et is the product of an elliptic and 
astigmatic Gaussian beam by an opaque disk of diameter D, then 

E T (£, n) = exp [c^ 2 + c n r, 2 ] ■ [1 - T(£ - £ 0)f7 - %)] • (2) 

V v ' " v ' 

—E(^.r]) shifted object 

The complex coefficients and are 

1 • TT _ _ 1 _ . TT 

For an opaque disk centered at the origin O, the transmittance function T(£, 77) in the object plane 
is: 



1, 0< y/e + v 1 < D/2, 
T (£> V) = < 1/2, < VF+F = D/2, ( 4 ) 



0, yji 2 + V 2 > D/2 > 0. 
^From Eqs. ([T]) and ([2"]l. the expression for A(x, y) can split into two integral terms, denoted A\ and 
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A 2 so that: 



A(x,y) = eMi ^ z) [A,-A 2 }. (5) 



The expression of these two terms is given in the Appendix (1A"1) . The complex amplitude A\ repre- 
sents the propagation of the incident beam without diffraction by the particle and A 2 contains the 
diffraction by a pinhole of diameter D. After some development detailed in the Appendix fA]), the 
expression of A 2 becomes: 

ttD 2 °° 
A 2 = -_cxp[$(£ ,%)] -cxp(ic z [(x ~ to) 2 + (y - Vo) 2 ]) • ^H) fe £ fc T fe (r, 7 ) cos(2fcfl), (6) 



Xz 

k=0 



with 



T fc (r, 7 ) = E^W^W^), (7) 

p=0 



and 

00 m+p 

^+2 P (r,7)=e X p(i7/2)E( 2m + 1 ) im ^(7/2)- E (-lj'^i^^™. (8) 

171= /-max(0,m-2fc-p,p-m) 

The coefficients w m / are given explicitly in [TT] and [T3j. Finally, the expression of A\ contains only 
the characteristics of the incident beam and A 2 contains a shifted linear chirp function linked to the 
object shift in the EAGB. The previous function is modulated by a series of Bessel functions which 
constitutes the envelope of the amplitude distribution of A 2 . 

A similar expression has been established in [11] for a particle located on the axis of the beam, 
i.e. £0 = and 770 = 0. This last expression is more general and can be applied to a particle located 
everywhere in the field of the beam. 

In our previous study, we demonstrated that a fractional-order Fourier transformation allows a 
particle located on the axis of the beam to be reconstructed. Unfortunately, the expression of the 
diffraction pattern (and particularly A 2 ) is so complex that it cannot be proved theoretically that a 
digital (or optical) reconstruction leads effectively to an opaque disk function. The demonstration is 
empirical. In addition, nothing proves that the fractional orders used for the reconstruction would 
be appropriate if the object were shifted transversely. 



In the next developments, we will thus prove theoretically that the reconstruction leads effectively 
to the object. We will then demonstrate that the reconstruction process is successful whatever the 
transversal position of the particle in the field of the beam. In addition, wc show that the same 
fractional orders have to be applied. We will thus be able to conclude that reconstruction is possible 
for a whole wide-field object or wide particle field in such a beam. 

Now, to obtain the desired accuracy, we should analyze the number of terms necessary in the 
series over k in Eq. ^ and over p in Eq. (J7J. The series which should be analyzed are in Eq.©. 
The first upper bound that is relevant here is [T1HT5] 



J i 



•(i) 



< 



i . A mva| 7 /4|' 

—7- mm 1,1 — 

(2m + 1) 1 / 2 V 



(9) 



The coefficients uj m i verify u m i > and u> m i = 1, so the w's are completely innocent. Using 
the bound © for the function Jm^) requires a rough estimate of the variable 7. In practical 
experiments, we have D « 10 _4 m, A « 10 _6 m, z « 10 _1 m, R q w 0.5 ■ 10 _1 then 7 sa 0.785 • 10 _1 + 
0.816 ■ I0~ 3 i. The right-hand side of Eq. ([9]) is less than 0.142 • 10 _1 for m > 1. The second upper 
bound we use concerns the Bessel functions J2&+2J+1- From [17], 9.1.62 on p. 362 and [M] we obtain: 

^2fe+2/+l(?') 



1 |r| 2fe+2 'exp(|r| 2 ) 
< min 1, — 



2 (2fc + 2Z + 1)! 

The evaluation of the accuracy concerns the product of the left-hand side of Eqs. (jSJ) and (|10[) 

•hk+2{m+p) + l{ r ) 



(10) 



d) 



(11) 



With the previous values, all the quantities in ([11)1 are less than 0.142 • 10 _1 for all (k,p,m) > 1. 
Thus, we only consider the case where fc = p = m = 0. In conclusion, the function V^{r, 7) is given 
with good accuracy by: 

F °(r, 7 ) ~ exp(i 7 /2)io(7/2)^. (12) 

r 

Note that this result is true if the object is far from the beam waist. The amplitude A 2 then becomes 



A2 = — /3 (<5)exp[$(e ,r, )]exp( lCz [{x - &) 2 + {v - Vo) 2 } ) • V °(r, 7). 



(13) 



2. A. Intensity distribution of the diffraction pattern 

The intensity distribution of the diffraction pattern in the quadratic sensor plane, denoted /, is 
evaluated from the Eqs. (TT5)) and (pi6"|) in the following way: 

I = AA = [A 1 -A 2 ][A 1 -M\ = [|Ai| 2 + |A 2 | 2 ] -2^{A{A 2 ), (14) 

where the overhead bar denotes the complex conjugate, 5ft denotes the real part. Thus, the intensity 
distribution recorded by the CCD sensor is described by the Eq. (JUJ). As one sees from the form of 
/, the first and second terms, i.e., |Ai| 2 and |^b| 2 , do not generate interference fringes with a linear 
instantaneous frequency, denoted fi(x), in the CCD plane, i.e. [18]: 

MX) = 2^— = dx = °" (15) 

But the third term exhibits a phase which is composed of a constant and a linear instantaneous 
frequency. This fact is important because the fractional Fourier transform is an effective operator 
for analyzing a signal containing a linear instantaneous frequency (linear chirp functions). From Eq. 
(Ti"4|) we write 

A\M = \A{A 2 \ exp [i arg (A X A$\ , (16) 
where arg = <j) — 4>q with 

<f> = c z (x 2 (Ms - 1) + y 2 (Al n - 1)) + 2c z (x£ + y Vo ) - arg f , (17) 

and 

0o = ^ + 3($(£o,r?o)) + arg (jo(|)) + arg ($(*)) - arg (K C K V ) , (18) 

where 3 represent the imaginary part of a complex number. The first term in (|17p yields a quadratic 
phase and the second term yields a linear phase. Remember that the aim of the reconstruction by 
means of FRFT is precisely to analyze a linear chirp. 

To give two different examples, it is necessary to fix the values of the parameters (wj, uj v ), (i?£, R^) 
and (D, A, z). 
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In the first case, the values are defined by (7mm, 1.75mm) for the beam waists, (—00, —50mm) 
for the wave's curvatures and the diameter D of the particle is equal to 150m and located at 
120mm from the CCD sensor. The wavelength of the laser beam is 632. 8nm. The distance between 
the cylindrical lens and the particle is 5 = 250mm. The particle is shifted from the origin by 
(£o)?7o) = (0.5mm, 0.2mm). Figure @ illustrates the diffraction pattern which is recorded by the 
camera. Note that the shift of the diffraction pattern (xo,yo) observed in the camera plane is not 
equal to the object shift (£oi??o)- If the particle is considered far from the beam waist then we have 
the formula: 

I A I ±z 

yo = -jxy- ( 19 ) 

The sign of z depends on the position of the particle compared with the position of the waist. If the 
particle is after the waist, the sign is positive. If it is in front of the waist, the sign is negative. Along 
cc-axis, the parameter A is infinite so that in the plane of the camera xq = £0 = 0.5mm and along 
y-axis, A = 50mm thus yo = 0.68mm. 

Now, in the second case, the values are defined by (7mm, 1.75mm) for the beam waists, 
(—00, 50mm) for the wave's curvatures and the diameter D of the particle is equal to 150m and 
located at z = 120mm from the CCD sensor. The distance between the cylindrical lens and the 
particle is S — 150mm. The particle is shifted from the origin by (£oj Vo) = (0.5mm, 0.2mm). Figure 
([3]) illustrates the diffraction pattern which is recorded by the camera. In the plane of the camera, 
a^o = £0 = 0.5mm, A = 50mm thus, from Eq. (fT9| . it leads to yo = —0.28mm. The diffraction 
pattern changes from elliptical fringes to hyperbolic fringes. These diffraction patterns will be used 
to reconstruct the image of the particle by FRFT. 
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3. Fractional Fourier transformation analysis of in-line holograms 

3. A. Two-dimensional Fractional Fourier transformation 

FRFT is an integral operator that has various application in signal and image processing. Its math- 
ematical definition is given in Ref. |19H21j . The two-dimensional fractional Fourier transformation 
of order a x for x-cross-section and a y for y-cross-section with < \a x \ < tt/2 and < \a y \ < ir/2, 
respectively, of a 2D- function I(x, y) is defined as (with a p = 

.^ ax ,a y [I(x,y)}(x a ,y a ) = / N ax (x,Xa)N ay (y,y a )I(x,y)dxdy, (20) 



where the kernel of the fractional operator is defined by 



N a Ax,x a ) = C(a p )exp «tt^- exp , (21) 

s^tana p y \ s^sinap ' 



and 



_ cxp(-i(fsigii(sinap)- ^)) 
(ttpj " l^sina p |V2 



Here p = x,y. Generally, the parameter s p is considered as a normalization constant. It can take 
any value. In our case, its value is defined from the experimental set-up according to |22j 

4 = N p -%. (23) 

This definition is presented in the Appendix [Cj N p is the number of samples along the x and y axes 
in both spatial I(x, y) and fractional domains. The constant 5 P is the sampling period along the two 
previous axes of the image. In our case the number of samples and the sampling period are the same 
along both axes, so that the parameters s p are equal to s. The energy-conservation law is ensured by 
the coefficient C(a p ) which is a function of the fractional order. One of the most important FRFT 
feature is its ability to transform a linear chirp into a Dirac impulse. Let g(x) = cxp(i7rxa; 2 ) the 
chirp function to be analyzed. If we consider — = — tana, and using the fact that: 

lim -^exp (~\ = 5(x), (24) 



£->o Wine V i 



then the FRFT of optimal order a of g(x) can be written as: 



& a [g(x)} =S(x a ). 



(25) 



Finally, by choosing the adequate value of the fractional order, a pure linear chirp function can be 
transformed into a delta Dirac distribution. 



3.B. Reconstruction: optimal fractional orders 

To reconstruct the image of the particle, the FRFT of the diffraction pattern (Eq. (|H)) ) must be 
calculated: 

#w«,[i] =-*wh [\Ai\ 2 ]-^, ay [2|^|cos(0-4>o)] +& ax , ay [\A 2 \ 2 ] (26) 

The terms |^4i| 2 and |^2| 2 do not contain a linear chirp, so they do not have any effect on the 
optimal fractional order to be determined. But the second term, denoted St contains a linear chirp. 
It will be considered for the image reconstruction of the particle. By noting that 2cos(0 — <po) = 
exp(—i(4> — <f>o)) + exp(i((j) — 4>o)), the second term of Eq. (f26)) becomes : 



&a x ,a y [2\AiA 2 \ COS(0 - 0o)l = CXp ( Z7T n~~ 1 ex P ( " 

" 1 J \ s 2 tana x / \ 



vl 



s 2 tan a y 



{/_+/+} (27) 



with 



I± = C{a x )C{a y ) / / \A 1 A 2 \cxv[i(cf> a ± 



)))] exp 



2m 



x a x y a y 
sin a x sin a. 



dxdy (28) 



The quadratic phase term of the FRFT is denoted by <j) a — (x 2 cot a x + y 2 cot a y ) . Let us recall 
that the FRFT allows us to analyze a linear chirp. Thus, if the fractional orders check the following 
conditions 



7T cot aZ pt 



c z (M x - 1), 



7rcota° p * 

= c z (M y - 1), 



(29) 



then, the FRFT of 7_ is no more than a classical Fourier transformation according to 



i- = x-J r 



Mr) 



■ exp 



(u,v), 



(30) 



with X = ^C(a x )C(a y )K i K 11 exp[<i>^ ,r]o)}l3°{S)exp(q)jo(^). The operator T is the 2D- 
Fourier transformation. The spatial frequencies u and v are equal to: 



CzVO 



s 2 sin(aS p ) 



sin(a2 pI ) * 



(31) 



As the Fourier transform of the product of two functions is equal to the convolution of their trans- 
forms then 



\Ji(r)~ 








r 





(32) 



Remember here that the variables p and r pertain to the same coordinates [x, y). With the shift 
theorem for the Fourier transform, the Hankel transform and the discontinuous Weber-Schafheitlin 
integral in [ [17]], 11.4.42 on p.487, we have: 



\Ji(r)~ 






r 





exp [— i2n(uX + vY )} x < 



1, 0< Vu 2 + « 2 < 



1/2, < Vu 2 + v 2 = ( 33 ) 



0, I > 0, 

with X = £ (1 — ic£/c z ) and Yq = rjo (1 — ic v /c z ). The function defined by the right-hand side of 
Eq. (|33[) , with spatial coordinates {x a ,y a ), has the aperture of the pinhole with diameter equal to 
the diameter D of the opaque particle. We have thus demonstrated our result: reconstruction with 
a FRFT leads exactly to the object. In addition, shifting the object does not modify the fractional 
order. This point is important because in the case of a particle field, a single fractional order couple 
along the x— axis and y— axis is necessary to reconstruct the particle image. If one wishes to 
determine the shift of the diffraction patterns in the (x a , y a )-plane, the coordinates that should be 
considered are: 

s 2 (u _ 2lk ) tan of*, s 2 {v - ^) tan<P^ (34) 

TT 7T y J 

This correction is necessary because the fractional Fourier transformation is not invariant by trans- 
lation: a space-shift in the spatial domain will lead to both frequency-shift and phase-shift in the 
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FRFT domain. The shift rule of the FRFT is expressed as: 

&a [f( x — b)](u) = cxp (inb 2 sin a cos a) exp (—i2nbu sin a) & a [f(x)](u — b cos a) (35) 

Note that the nature of the Gaussian beam implies that the reconstructed image of the object is 
convolved by a 2D Gaussian function. The sign of the object function (opaque object with a luminous 
background) is retrieved by applying an inversion of 7_ that is realized by the minus in front of the 
second term of Eq. (|26l) . 

3.C. Numerical experiments 

The simulations of the particle image reconstruction are realized from the diffraction patterns illus- 
trated by the Figs. ([2]) and (|3|). The diffraction patterns consist of a 512 x 512 array of 11m x 11?ti 
size pixels. Consider the diffraction pattern presented in Fig. ^ and produced by a particle of 
D = 150m in diameter located at z = 120mm. The optimal fractional orders obtained from Eq. (|29|) 
are a° pt = —0.564 and a° pt = —0.850. The image of the reconstructed image is shown in Fig. In 
this representation, the squared modulus of the FRFT, i.e. \J^a x ,a y [I] \ 2 , is taken. The shape of the 
particle image is not modified: the width of the Gaussian function in the Eq. (|32|) is greater than the 
diameter D of the particle (typically 679m along x-cross axis and 490m along y-cross axis). Now, the 
reconstruction of the particle image from the diffraction pattern illustrated by Fig. ([3]) , is realized 
by a fractional Fourier transformation of optimal orders a° pt = —0.564 and a° pt = 0.664. Figure 
(J3J) illustrates the result of the reconstruction. In both cases (Fig. (0| and Fig.©) reconstruction by 
FRFT is successful. Note that we have checked that the apertures resulting from Eq. (|3"2"|) and the 
image of the reconstructed particle give the same diameter D. 

3.D. Experimental results 

As the reconstruction process is successful whatever the transversal position of the particle in the 
field of the beam, for the same fractional orders we can now carry out a reconstruction for a wide-field 
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object. 

The previous theoretical developments and numerical experiments have been tested by using an 
RS-3 standard reticle (Malvern Equipment). This reticle is an optical glass plate with a pattern of the 
word "ELECTRO" photographically deposited on the surface. The word "ELECTRO" spread over 
6.5mm. The reticle is located at S rj 156mm from the cylindrical lens (CL). The distance z between 
the CCD camera and the reticle is approximatively equal to 117mm. The opaque word "ELECTRO" 
is in front of the waist of the beam. The intensity distribution of the EAGB diffracted by the word 
"ELECTRO" is shown in Fig. ©. The image of the word "ELECTRO" is reconstructed by the 
fractional Fourier transformation. Fractional orders have been adjusted to obtain the best contrast 
between the reconstructed "ELECTRO" word and the background. Thus doing, the approximate 
orders are a° pt = 0.505 and a° pt = —0.785. The image of Fig. ([7]) shows that two previous orders 
allow all parts of the image of the object to be reconstructed. Note that the word is reversed along 
the y— axis. This is due to the Gouy phase shift of a Gaussian beam that propagates from — oo to 
+oo through the focus point. In this case, the phase shift predicted is tt. 

4. Conclusion 

The effect of a transversal object shift in an elliptical, astigmatic Gaussian beam does not affect 
the optimal fractional orders required to reconstruct the image of a particle or any other opaque 
object in the field. In this publication an analytical model has been developed to prove this. In 
this development, an opaque particle is considered. Experimental result presented by the recording 
and reconstruction of the word "ELECTRO" show that our method works, even when the object is 
spread out over the whole image field. 
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A. Appendix A : Expression of the terms A\ and A2 

A. A. Expression for the amplitude distribution A± 

The development of integral A\ has been given in a previous paper by [TT] 

A 1 = K^K n exp (-^- P T Np ) j exp (i^p T Mp) 



(36) 



where p T represents the vector (x y) and the factors K q with q = £, r\ in Eq. (|36[) are dchncd by 



K„ 



and the diagonal matrices N and M by 



Xz 



1/2 



(37) 





( N x ^ 




( M x ^ 










N = 




M = 






^ Ny J 




y M y ) 



(38) 



with 



N q = n 



Xz 



1 + 7T S 



1 



M q = 1 + 7T 



1 



1 + - 2 T^ t- 1 



(39) 



v4.Z?. Expression for the amplitude distribution A2 

To develop the second integral of A2, involving the product of E(£, rf) and T(£ — £0, r\ — rjo), i.e.: 



A2 = exp[^+^) i 



we first replace £ by £ + £0 an d rj by + 770 and get: 



T (£ "b V J 
Xz 



exp 



-1— [xi + yn) 
Xz 



d£dr], 
(40) 



exp \f-(x 2 + y 2 )] f 
M = \ z J E ^ + ^,V + Vo) exp 



f- ((e + £o) 2 + (^ + r,o) 2 ) 



exp 



27T 



d£,dn. (41) 



The domain I? must be defined as the disk with center and diameter D. By considering that 
c z = tt/(Xz) and by restating A2 in cylindrical coordinates as follows: £ = Da cos(</?)/2 and r/ = 



13 



Da sm(ip)/2 for the object plane, we obtain: 



D 2 

A2 = JX~ Z ex P [ c ^o + c vVi + ic z [(x - Co) 2 + {y~ Vo) 2 ]] 



I I exp [«7cr 2 ] exp \i5a 2 cos(2</?)] exp [iacr cos ip + ibasimp] adadip (42) 
Jo Jo 



with 



D 2 .D 2 . . . ,D 2 . 

7= — c 2 -«— (c c 4-c,,), 6 = i — (Cr, - c 6 ), 



a = Dc z [Co (1 - «c 5 /c z ) -x], b = Dc z [rj (1 - ic n /c z ) - y] 



(43) 



By writing: 

a cos <p + b sin (p = r cos(tp — 6) (44) 

for which we have the condition 

a = r cos 6, b = r sin 6 (45) 

with complex r and 6. This representation is discussed in some detail in Appendix [B] By means of 
the following equalities in [ [17|], 9.1.41 - 9.1.45: 

+ 00 

exp [iSa 2 cos (2<p + 26)] = J (5a 2 ) + 2^i k J k (5a 2 ) cos 2k((p + 6), (46) 

fe=i 



and 

1 r 2 ^ 

exp(in9) exp[ix cos 9}dd = i n J n (x), (47) 



n 



2- 

the expression of A 2 becomes: 

ttD 2 °° 
A 2 =1L— exp [$(&,%)]■ exp (ic x [(x - Co) 2 + {y - Vo) 2 ]) ■ T k {r, 7) cos(2fc#), (48) 

fc=0 

with £fc = 1/2 if fc = and 1 otherwise. The parameter denoted 3>(£o)77o) is equal to [c^Co + c i]Vo] ■ 
The function T k (r, 2 7 ) is defined as: 

oo 

T k (r, 7) = E $k + 2p( S )V 2 f +2p (r, 7 ), (49) 

p=0 
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where the coefficients /3f£ +2p are given by the analytical development of Tf. in Appendix of [ [TT|] . 
Recall here that the expression of V 2 2 £ +2p (r, 7) is: 

00 m+p , . 

V^ +2p (r,7) = exp(i 7 /2) ^ (2m + l)i m j m ( 7 /2) ■ £ (-l)W J2fc +^ +1 ^ . (50) 

m=0 Z — max (0,m — 2k—p, p — m) 

The coefficients uj m i is given explicitly in [ [TT]] and [ [Tl] ]. 

B. Appendix B : Elaboration of condition (|45p 

With u = a + ib, v = a — ib, we should find r and to = exp(i0) such that 

TLU = U, t/lU = V (51) 

Assume that ti ^ 0, w ^ 0, and write u = r expect), v = s exp(z/3) with r, s > and a, /3 € R then 
we see that uj = exp(i#) given by 

1/2 

t = (rs) 1/2 exp(i(a + /?)/2) , w = exp(i0) = exp (i(a - /3)/2) (52) 

satisfy the relations (|51l) . This does not work in the case that u = or v = 0. Indeed, when a = 1, 
b = i we get from (|45[) and cos 2 9 + sin 2 6* = 1 that r = a + ib = 0, i.e., r = 0. 

C. Appendix C : Definition of s p 

To determine the value of s p , it is necessary to write the definition of the one-dimensional fractional 
Fourier transformation in the particular case of a = ir/2: 

T 7r/2 [I(x))(x a )=C(ir/2) I(x)exp(~i2ir^) dx. (53) 

J-00 ^ s / 

Its discrete version is 

^ /2 [I(m)}(k)=C(n/2) J2 I(rn)exp(-i2ir m6x J 6x A s x , (54) 

m=-N/2 ^ ' 

where S x and S Xa are the sampling periods of I(x) and its transform. The sampling periods are equal 
to S x , and N is the number of samples of I(x) and its transform. The relation (|54j) can be written 
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as the discrete Fourier transformation of I(m) according to: 



N/2-1 

mk 



^ /2 [I(m)](k) = C^/2) J2 I(m)exp[-i2^)5 x . (55) 

m=-N/2 



By identification of Eqs. (|55|) and (fMl) . one obtains: 

so s 2 = NS x S Xa = N6'£. (56) 



SxSx ° 1 „2 _ ,\rc r _ A ,r2 



In the case of two-dimensional function one finally has s 2 , = N p ■ <5 2 with p = £, ?y. 
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Fig. 1. Numerical and experimental optical set-up. 
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Fig. 2. Diffraction pattern with , = 7mm, lo v = 1.75mm, = oo, R v = —50mm, 
D = 150m, A = 632. 8nm, z = 120mm, 5 = 250mm, £o = 0.5mm and rjo = 0.2mm. 
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Fig. 3. Diffraction pattern with , uj^ — 7mm, lo v = 1.75mm, = oo. R v = —50mm, 
D = 150m, A = 632.8?im, z = 120mm, 5 = 150mm, £o = 0.5mm and rjo = 0.2mm. 
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Fig. 7. Fractional Fourier transform of the diffraction pattern with a° pt = 0.505 and a° pt 
-0.785. 
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